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We study the depinning transition of a driven chain-like system in the presence of frustration and 
quenched disorder. The analysis is motivated by recent transport experiments on artificial vortex- 
flow channels in superconducting thin films. We start with a London description of the vortices and 
then map the problem onto a generalized Frenkel-Kontorova model and its continuous equivalent, 
the sine-Gordon model. In the absence of disorder, frustration reduces the depinning threshold in 
the commensurate phase, which nearly vanishes in the incommensurate regime. Depinning of the 
driven frustrated chain occurs via unstable configurations that are localized at the boundaries of 
the sample and evolve into topological defects which move freely through the entire sample. In 
the presence of disorder, topological defects can also be generated in the bulk. Further, disorder 
leads to pinning of topological defects. We find that weak disorder effectively reduces the depinning 
threshold in the commensurate phase, but increases the threshold in the incommensurate phase. 



PACS numbers: 74.60. Ge, 71.45.Lr, 83.50.Lh 
I. INTRODUCTION 

The depinning transition in driven systems has at- 
tracted a great deal of attention recently. The phe- 
nomenon can be observed in charge, density wave sys- 
tems submitted to an electric field,EJ in magnetic butt, 
bles moving under an applied magnetic field |#ssadient,EI 
in current driven Josephson junction arrays,tru and in 
the two-dimensional electron gas in a magnetic field, 
which at low densities forms a,Wigner crystal, but can 
move under an applied voltage.Q Depinning is also impor- 
tant to understand tribology and solid friction,El surface 
growth of crystals with quenched bulk or substrate dis- 
order, and £he dynamics of domain walls in incommensu- 
rate solids. u A very prominent example of driven systems 
displaying a depinning transition are vortices in type-II 
superconductors.El 

The H — T phase diagram of type-II supercondwkj 
tors displays regions of both elastic and plastic AowBlI 
The interplay between the elasticity of the vortex lat- 
tice and the impurities present in the substrate leads 
to a rich phenomenology with many static and dynamic 
phases. A crucial question in both the dynamics and 
statics is whether - in addition to thermal fluctuations 
- quenched disorder produces topological defects in the 
periodic structure. Whereas in the absence of topolog- 
ical defects it is sufficient to consider only elastic de- 
formations with depinning causing elastic flow, the dy- 
namics will be governed by plastic flow, if topological 
defects exist. One expects plastic motion to become im- 
portant for either strong disorder, high temperature, or 
near the depinning transition in low dimensions. In these 
cases, depinning is observed to occur through "plastic 
channels" between pinned regions. This type of plastic 
flow has been found in numerical simulations of a two- 
dimensional thin film geometry.liHtL3 Above the thresh- 



old, the filamentary channels become both denser and 
broader when the driving force is increased. Measure- 
ments of the differential resistance of MoGe films dis- 
play abrupt steps, which could be interpreted in terms 
of plastic depinning. t3 Other experimental observations 
that haye . been attributed to plastic-,flow are the peak 
effectjLaEJ vmsual broadband noise,ll3 and steps in the 
I-V curve.EZI 

In order to study the plastic depinning of vorticeSy-afc. 
tificial easy-flow channels have been manufactured£3'E£l 
The samples are typically made of type-II superconduct- 
ing thin films, which consist of a fairly strong pinning 
layer on top of a weakly pinning base layer. The artifi- 
cial vortex flow channels are fabricated by etching away 
stripes of the stronger pinning top layer. In the Shub- 
nikov phase, vortices penetrate the sample. Applying 
a current perpendicular to the channel direction, the re- 
sulting Lorentz force in the direction of the channel drives 
the vortices. Since the influence of point-like material 
defects in the weak pinning channel is negligible, chan- 
nel vortices are mainly pinned indirectly via the inter- 
action with the stronger pinned vortices in the channel 
environment. Above a threshold force, plastic depinning 
of the vortices inside the channel takes place. In con- 
trast to natural channels, the depinning threshold force 
displays interesting oscillations when the externally ap- 
plied magnetic field is varied. The position of the thresh- 
old force maxima hint at commcnsurability effects be- 
tween the. vortex lattice in the environment and in the 
channel.LJ The magnitude of the depinning force minima 
and maxima indicate that lattice distortions produced by 
quenched disorder in the pinned channel environment is 
relevant in these samples.EJ 

In this work, we present theoretical studies of the de- 
pinning and dynamical behavior of vortices in artificial 
easy-flow channels. Frustration and disorder are essen- 
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tial ingredients for determining the depinning transition 
in this driven system. Motivated by recent experiments 
on vortex-flow channels, we develop a model which takes 
into account inhomogeneities, such as the sample bound- 
aries and disorder. We find that depinning occurs due to 
boundary or bulk nucleation, respectively, in the absence 
or presence of disorder. Further, we clarify the role of 
disorder and frustration and compare our results to the 
available experimental data. 

First, we study this problem considering a perfectly 
ordered vortex-lattice in the environment. Starting 
from a London description of vortices inside and out- 
side the channel, we derive the-,coefficients of a general- 
ized Frenkel-Kontorova model.E3 The force- velocity curve 
of this simple model displays a drastic decrease of the 
threshold pinning force when topological defects enter 
the system via the system boundaries. However, the ex- 
perimental data show a smooth variation of the pinning 
force as a function of the magnetic field. The behav- 
ior suggests that disorder in the channel edges may lead 
to smoothening of the transition from the topologically 
ordered to the defective state. In order to investigate 
the transition in presence of disorder, we take into ac- 
count that the vortex lattice in the channel environment 
is distorted due to impurities in the superconducting sub- 
strate. To determine the characteristic current-voltage or 
force-velocity curves, the problem is most conveniently 
tackled with numerical molecular dynamics simulations. 
This method enables us also to nicely visualize how de- 
pinning occurs in the topologically ordered and defective 
regimes. It turns out that systems that are topologically 
ordered in the static phase depin at the boundary or by 
generating mobile defect droplets at weak spots in the 
random channel potential. In the defective static phase, 
the pre-existing topological defects are pinned by the ran- 
dom potential. Then, the depinning transition occurs via 
a release of these defects at the threshold. 

The paper is organized as follows: In Section II, we 
discuss the experiments that motivated the theoretical 
studies; in Section III we introduce the model and present 
the results concerning the depinning transition. Finally, 
the discussion and conclusions are presented in Section 
IV. 



II. EXPERIMENTAL MOTIVATION 

Let us start with a short overview of the experimental 
status of art. In order to obtain a more detailed under- 
standing of the depinning process in disordered supercon- 
ducting films, .experimentalists fabricated artificial vor- 
tex channels. A sketch of a typical device is presented 
in Fig. |l|. The devices are manufactured by first absorb- 
ing a thin layer of NbN (strong pinning material) on top 
of a NbjGe film (weak pinning) and then etching away 
a few hundred channels of width w from the top layer. 
In the strongly pinned region, the vortex-lattice spac- 
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FIG. 1. Schematic plot of the artificial flow channel ge- 
ometry. The channels are manufactured by etching away the 
strongly pinning top layer material (grey area) in stripes of 
width w. Applying a magnetic field H generates vortices. In 
the strong pinning area the vortices are static. Vortices in the 
weakly pinning channel are mainly pinned due to the interac- 
tion with the strongly pinned ones in the channel edges (grey 
area). In presence of a current density j a Lorentz force fz, 
acts on the vortices. At the depinning transition the channel 
vortices start to flow whereas the strongly pinned vortices in 
the channel environment remain static. 



ing is a = (4/3) 1 / 4 ($ /S) 1 / 2 , such that a b = <f> /B, 
with $o the flux quantum and bo = (-\/3/2)ao, see Fig. 
2a. In order to measure the effect of a single vortex row 
in the channel, the channel width w is chosen to be of 
the order of ao. Therefore, the vortex-lattice spacing in 
the channel is a = <& /wB. Depending on the mismatch 
parameter \ — (o/ao) — 1 = (bo/w) — 1, a commensurate- 
incommensurate transition can be observed. In the ex- 
periments described here, ao ~ w ~ 100 nm and the 
distance between the channels is 10 /xm. Due to the ma- 
terial steps at the channel edges, screening currents flow 
along the surface which interact with the vortices in the 
channels. As a consequence, we can expect the first rows 
of strongly pinned vortices near the channel edges to be 
lined up at fairly straight along the channel. Therefore, 
we can assume that the transversal shifts are negligible. 
On the other hand, longitudinal displacements may occur 
due to the presence of impurities in the superconducting 
material [see Fig. 0(b)] . Here we will concentrate on the 
longitudinal dynamics of the vortices within the channel. 

In the channel, pinning due to impurities is very weak. 
Below the depinning threshold, vortices in the channel 
are localized due to the interaction with the vortices in 
the static environment. In the presence of a current den- 
sity j, a Lorentz force fx, acts on the vortices. As long 
as the vortices are pinned, the supercurrent is flowing 
dissipation-less through the sample. At a threshold force 
f c the channel vortices depin and start to move. The 
vortex motion causes dissipation of energy and a voltage 
drop occurs across the sample boundaries. Hence, a well- 
defined critical current density j c can be determined by 
measuring the current- voltage characteristics of the sam- 
ple. Fig. shows shear force density data deduced from 
current-voltage measurements.Ea The oscillations in the 
shear force density F s as a function of B are related to 
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FIG. 2. (a) Easy flow channel environment (shaded) in the 
absence of disorder. Vortex positions in the environment are 
marked by open circles, (b) Immobile vortices (dots) in the 
environment in the presence of weak disorder and one row of 
mobile vortices in the flow channel. 
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FIG. 3. The experimentally obtained critical shear force 
density F 3 versus magnetic field B for a channel sample with 
w « 290 nm at a temperature T = 1.7 fM 



the (in)commensurability of the vortex lattice &o with the 
channel width w. Naively, we would expect the maxima 
of F s to occur at integer values of the ratio w/bo, where 
the vortex lattice spacing inside the channel is commen- 
surate with the width of the channel. F s should then be 
reduced for non-integer w/bo. This could qualitatively 
be explained by the development of misfit dislocations 
along the channel edges Jialij Although this intuitive pic- 
ture holds for the case of weak (or zero) disorder, these 
features change drastically in the limit of a strongly dis- 
ordered environment. Indeed, as it will be explained in 
Ref. |Ea], in the latter limit, the minima and not the max- 
ima of F s correspond to an integer number of vortex-rows 
in the channel. In this paper, however, we restrict our- 
selves to the limit of weakly disorderd environment and 
investigate the competition between boundary and bulk 
nucleation for activating the depinning process. 



III. A MODEL FOR ARTIFICIAL VORTEX 
CHANNELS 



model for vortices in an artificial easy-flow channel. 



A. Channel with a perfect vortex lattice in the 
environment 

Let us first consider a row of straight vortices that are 
aligned along the a;— axis at positions maoe x , where m 
is an integer, ao the distance between adjacent vortices 
and the unit vector in the x— direction. In a type- 
II superconductor within London theory the interaction 
energy between two straight vortices of length I at r and 
is U(t) = U K (r/X), where U = $gZ / (8ir 2 A 2 ) . Here A is 
the penetration depth, $o = hc/2e is the flux quantum, 
and I ^> A the sample thickness. Then the potential 
energy E°(r) that is felt by a test vortex at position 
r = (x, y) interacting with a row of vortices placed at 
nboe y is 



U(r - R m ,„) 



(1) 



where we introduced the lattice vectors R m , rl = ([m + 
n/2]ao,nbo) of a hexagonal lattice. Fourier transforming 
and using the Poisson sum rule the latter can be recast 
into 



cos[g„(a; — nao/2)] 



2na 



o 



dk v Ue ik y (v - nbo) 



(2) 



where q v = 2iri//ao with integer v and U is the Fourier 
transform of the potential, which is given by 



U(q v , k y ) = — 



2irU 



kl + A- 2 ■ 



(3) 



Performing the integral over k y leads to 

E° n {r) =Y,M\V - n&ol) cos[q v {x - na /2)}, (4) 



where 



Bu(y) 



ttUq 
aoql 



(5) 



and q' v = y/qJ+X 2 . 

We construct an easy flow channel by building a two 
dimensional vortex lattice, but leaving a region of width 
w along the x— axis empty (see Fig. ga). 

If we consider a hexagonal vortex lattice in the channel 
environment, the potential in the channel is 



E oc (r) = J2 [K (r - b/2) + E°_ n (r + b/2)] 



(6) 



Starting from a London description of vortices in a 
type-II superconductor, we now derive a simplified Id 



where b = (w — bo)e y . Summing over the vortex rows n 
one finds 
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Kc(r) 



]A v (y) cos (q„x) 



the Fourier coefficients A u (y) read 

2cosh(^y)S v (3£=^ 



(_l)v e gtbo _ 1 



(7) 



(8) 



Since the magnetic inductions B used in the experi- 
ments cover the entire range up to the upper critical one 
B C 2, a more general expression for A v than the London 
limit that has been discussed until now is needed. How- 
ever, the theory presented here can easily be extended 
to larger magnetic inductions. First one takes into ac- 
count the finite diameter of the vortex core which is of 
the order £ in the vortex-vortex interaction potential, 
U — > Uo[Ko(r/X) — Ko(r/£)]. Second, one replaces A — » 
A' = A/(l - B/B c2 f/ 2 and = £/(2 - 2B/B c2 f/ 2 

to take into account the reduction of the superconduct- 
ing order parameter at large magnetic fields £-3 We then 
obtain 



Au(y) = 



ttU (1 - B/B c2 ) 



«o 



2cosh(g^)e-«^( u '- bo )/ 2 
ql [(-l)"e^o - 1] 



2cosh(^'j/)e-«" ! '( u '- bo )/ 2 
q'l [(-l)"e9" b o - 1] 



(9) 



where now ql = y/q 2 .+ (A') -2 and q" = \/g 2 + (£')" 



B. Equation of motion 

The overdamped dynamics of a vortex with index m 
at position r m inside the channel is described by 



■qr n 



f ■ e x - V£ oc (r„ 



J] VZ7(r m - r„), (10) 



where / = j$>o/c is the magnitude of the Lorentz force 
that drives the vortices in presence of a current density 
j. The viscous drag coefficient r\ is related to the flux 
flow resistivity pjf by 77 = B$o/[c 2 Pff(B)]. The sum is 
taken over the positions n of all other vortices inside the 
channel. 

In the simplest case the channel width is w ~ 60 
such that only a single vortex row is inside the chan- 
nel. Then the motion of the mobile vortices in the 
y— direction is essentially guided by the channel poten- 
tial, whereas the interaction between vortices in the chan- 
nel does not contribute significantly to the motion in the 
y— direction, r\y m ~ ~d v E oc {v). However, the motion 
in the x-direction is determined by both the interactions 
between mobile vortices and the gradient of the chan- 
nel potential. To simplify matters, we neglect the mo- 
tion in the y— direction, y m — 0, such that the equa- 
tion of motion can be simplified to a one-dimensional 
one. Further, we restrict our considerations to A > do 



as in a typical experiment. Then, the amplitudes A u (0) 
fall off exponentially fast, ^4^(0) ~ Uo exp(— 7rv3i/), and 
the approximation to consider only the first harmonic 
q = qi = 2-k / a$ of the channel potential is a good one. In- 
troducing /i = 2g| j4i (0) I , and restricting the interaction 
between vortices in the channel to next-neighbor spring- 
like forces, the equation for the overdamped longitudinal 
motion reads 



r\x n 



dV 

QXrn. 



(11) 



where the potential energy of the vortices in the channel 
is given by 



V 



1 - [1 - cos (qx m )] - fx n 



T / 4 2 \ x m+n 



naY 



(12) 



It hasp- the form of a generalized Frenkel-Kontorova 
modelB3 The interactions between vortices inside the 
channel are approximated by Hookian springs with spring 
constants n n = (U /X 2 )K^(na/X), where the double 
prime denotes the second derivative. Notice that the dis- 
placements x are measured with respect to the vortex- 
lattice preferred position a in the channel. 

The Frenkel-Kontorova model has been intensively 
studied close to equilibrium, / ~ 0, see Ref. [Ell. The 
reduced dimensionless elasticity 



g= — > n K n 



e nV3 X 

»1 

07T a 



(13) 



together with the winding number £1 = a/ao crucially 
determine the behavior of the system. For rational Q the 
vortex chain is commensurate with the periodic channel 
potential whereas for irrational O it is in an incommen- 
surate state. The commensurate-incommensurate tran- 
sition is a continuous transition that occurs at finite mis- 
matches x (recall that \ — (°. — ®o)/ a o)i since the creation 
of discommensurations costs energy. 

If the discreteness of the chain is relevant, Peierls- 
Nabarro barriers may exist. The Peierls-Nabarro is the 
energy barrier that has to be overcome for a translation, 
x m — > a: m -|-i. Whereas this barrier is always finite for ra- 
tional Q , it may vanish in the incommensurate state: if g 
is lower than a critical value, g c (a/ao), the incommensu- 
rate state is pinned, however, for g > g c (a/ao), the Aubry 
transition to a sliding state with truly vanishing critical 
force takes place. c£l For pinned defective configurations 
the Peierls-Nabarro barrier, which determines the corre- 
sponding pinning force Jpn, depends on g. Large g ^> 1 
imply that an isolated defect having a size ~ Id extends 
over several lattice constants. Then, the Peierls-Nabarro 
barrier is nearly vanishingc3 and the pinning force is 



PN 



64ir 2 gu exp(-7r 2 v / g). 



(14) 
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Since fpN <C 10 _4 /x we neglect it in the following. 

Since we are most interested in the regime g 3> 1, 
it is convenient to study the model in the continuum 
limit. Further, we take into account the finite length 
of the system and consider open boundary conditions. 
Rewriting Eq. ( p^ ) in terms of the displacements of 
the vortices from the environmental lattice positions, 
u rn = x m — mao and substituting mao — > x, u m — > u(x), 
(u rn+n - u m )/(na ) -> d x u{x), J2 m ~* / dx/a Q , we ob- 
tain the continuum equation of motion 



rju(x) = —0,0 



SV[i 



Su(x) 
+ f + Xk[8(x 



_ d 2 u(x) 

dx 2 
L) 5{x)\ 



/zsin [(711(2:)] 



(15) 



where we have introduced k = 2iTao^g and the last term 
is an effective surface force that arises at the open sample 
boundaries in the presence of frustration (L is the channel 
length). The energy functional reads 



where the integral runs from to L. It can be decom- 
posed into 



V = V SG + V X + V Q , 



(16) 



where 



Vsu{<,}= I ^ jf (g) + £[l-cos(gu)]-/u 



is the energy functional of the sine-Gordon model, 



V x [u] = -—HL)-u(0)] 
a 



(17) 



is the frustration energy due to the mismatch \ deter- 
mined by the values of the displacement field u(0) and 
u{L) at the boundaries, and Vq is an irrelevant offset that 
is omitted in the following. 



C. Commensurate-Incommensurate Transition 

We now present a short review of the commensurate- 
incommensurate transition. Then, we will extend the 
picture to discuss the role of edge barriers for defects in 
finite systems, which is crucial to understand how dis- 
commensurations penetrate a sample in the absence of 
thermal or quantum fluctuations. 

The extrema of V[u] are found solving the variational 
problem SV/Su = 0. In the absence of frustration, \ — 0, 
and for / = 0, uniform static solutions exist. They are 
the stable, u Si „ = aon, unstable, u un = (2n + l)o /2, 
and kinked 



Uk,n(x; x c ) = u S!n + - arctan 



exp 



(18) 



solutions of the sine-Gordon model at zero driving force. 
Here, Id = aso-v/5 is the width of the kink. The 
corresponding anti-kink solution reads u a , n {x] x c ) — 
Uk,n{—x]—x c ). The commensurate-incommensurate 
transition is a transition that occurs in equilibrium (/ = 
0) when it becomes energetically favorable to have a fi- 
nite density of discommensurations in the system. Ne- 
glecting effects of system boundaries, the mismatch \d 
at which the commensurate-incommensurate transition 
takes place is then quickly found by comparing the energy 
of the kinked solution with the energy of the stable homo- 
geneous one V[ttfc >n ]-V[u a ,„] = (4/7t)(k/V<7) 1/2 -kx = 0, 



4 ( n 
Xa = - — 



7T 



1/2 



W9 



(19) 



Though at the commensurate-incommensurate transition 
Uk, n or u a ^ n have the same energy as u StTl , a barrier has 
to be overcome in order to make the transition from 
u s , n — > v-k,n or u s , n —> u a ,n, see upper curve in Fig. 4. In 
the absence of fluctuations, as considered here, discom- 
mensurations can only enter the system when this bar- 
rier vanishes (lower curve in Fig. 4) . In general, the bar- 
rier vanishes at a sufficiently strong frustration or driving 
force. We therefore define a threshold frustration xdj), 
which is a function of the driving force /. Finite driving 
forces / > are considered in Section [II D, where de- 
pinning is studied. In the following we investigate how 
(anti-)kinks enter the system at equilibrium / = and 
determine the zero-force threshold x c (0). 

Let us first determine the energy that is needed to de- 
form a uniform state into a kinked one. We note that due 
to frustration kinks can only be created spontaneously at 
the system boundaries, x = or x = L. Of course, in the 
presence of thermal fluctuations, quantum fluctuations, 
or quenched disorder, deep in the bulk kinks can in prin- 
ciple be created in the form kink-anti-kink pairs. How- 
ever, in the absence of fluctuations as considered here 
this is not possible. The reason is that kink-anti-kink 
pairs cannot gain frustration energy for a spontaneous 
kink emergence, since u(0) = u(L) — u s-n and hence 
V x = 0. Thus, (anti-)kinks can only enter the system at 
x = and x — L. In the following we discuss the pen- 
etration of an anti-kink at x = 0, having in mind that 
the same holds for x = L and for kinks. Note that in a 
finite system with boundaries at x — and x = L one 
has lim 2 ; e ^_ 00 u o> „(0; x c ) — u s>n . Hence, a deformation 
from Ws,n to iia.n can be achieved by pushing an anti-kink 
centered at x c from x c — — oo to x c > 0. The energy of 
the anti-kink solution relative to the uniform one as a 
function of center coordinate x c is then 

IW{x c ) = V[u a . n (x; x c )} - V[u a>n (x)] 

= a^R {x[w a: „(0; x c ) - u a , n (L; x c )] 
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FIG. 4. Potential energy of an (anti-)kink AV in units 
of Vk = kxd/2 as a function of the (anti-)kink center 
x c . Displayed are values close to the system boundary at 
x — for frustration parameters x! = Xa (upper curve) 
and |x| = Xc(0) (lower curve). At the classical commensu- 
rate-incommensurate transition, where x! — Xa, an entry 
barrier has to be overcome with the help of fluctuations to 
make a transition from u s ,n to Uk, n or u a , n - The entry bar- 
rier vanishes at xj = Xc(O) where (anti-)kinks can penetrate 
the system. 



of a boundary. Instead, the stable solution is a virtual 
anti-kink with a center x c localized outside the system at 
x Cj ± < 0. At the boundary, the chain thus tries to adapt 
optimally to the frustration to reduce its energy. For 
X = — Xc(0) one finds x c = 0, which means that half of 
the kink is already inside the system and that it can gain 
more energy by sliding towards the center of the system. 
The scenario is the same for an anti-kink entering the 
system at x = L. For a kink the description given above 
is identical except that \ > 0. The kink entry barrier 
vanishes at \ — Xc(0). 

The picture of the commensurate-incommensurate 
transition is thus drastically modified in the presence 
of system boundaries when there are no physical mech- 
anisms like thermal or quantum fluctuations that are 
needed to cross the edge-barrier. In fact, since in the ab- 
sence of fluctuations a system with boundaries remains 
commensurate for |x| < Xc(0), we identify the threshold 
a t Ixl = Xc(0) with the commensurate-incommensurate 
transition of a finite, purely mechanical system. 



dx 



dx 



(20) 



Since the tail of the kink falls off exponentially, we can 
neglect the influence of the boundary at x — L Id- 
For simplicity, we then consider a semi-infinite system, 
L — > oo, and obtain 



AV(x c ) 



«Xc 



4 x 



arctan 



exp 



1 



2 

+ "» h (S)}- < 21 > 

Minimizing AV with respect to x c we find a minimum at 

TTXc 



x cA = l d \n -■ 



and a maximum at 



2X 



4x 2 



1 



(22) 



x c ,2 = /din - 



At the frustration 
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X = -Xc(0) = -~Xa = p, 



(23) 



(24) 



the minimum and the maximum merge into a saddle 
point at x c = where the entry barrier vanishes and 
an anti-kink flows freely into the system. It is interest- 
ing to note that for x < the minimum of the anti-kink 
energy relative to the uniform solution is always nega- 
tive, AV(x c ,i) < 0. This means that in a frustrated sys- 
tem the uniform solution u s:n is unstable in the presence 



D. Depinning in the presence of boundaries 

In the following we investigate how the chain inside the 
channel actually depins in presence of a driving force, 
/ > 0. In the simplest case, for w = &o, we have a 
commensurate state without frustration, a — ao and x = 
0. The chain locks perfectly to the potential and the 
threshold depinning force is 



/c(0) = li. 



(25) 



If the system is frustrated, bo ^ w, depinning occurs 
via mobile discommensurations which in the sine-Gordon 
model are represented by kinks or anti-kinks. As in the 
equilibrium case, for finite driving forces discommensura- 
tions enter the system when their entry barrier vanishes. 
Whether they are mobile or not depends on further bar- 
riers that may exist in the bulk. 

In the continuous limit as discussed here, the defec- 
tive state is not pinned. Hence, in the continuum model 
the threshold Xc(0) indicates the change from a static 
equilibrium ground state (u s , n ) to a mobile one (v,k, n or 

Note that an exit barrier exists for a single kink at the 
other boundary (imagine the mirror image of the entry 
barrier as shown in Fig. |] at the other end of the sys- 
tem). However, the exit barrier becomes irrelevant in 
the presence of further kinks. This can be easily under- 
stood by the following argument. Suppose a kink enters 
the system, freely flows to the other end, and then be- 
comes trapped by the exit barrier. Then a second kink 
follows and interacts with the first one. If it would move 
"adiabatically" it would become trapped by the inter- 
action with the first kink, which mediates the pinning 
force. However, the first kink would experience the inter- 
action of the second kink, too. The resulting interaction 



G 



force is of the same magnitude as the pinning force, but 
of opposite sign. Hence, the total force is zero and the 
first kink is released. The second becomes pinned for a 
while until it is released by the third and so forth. For 
non-adiabatically moving kinks, the successor does not 
even become pinned by the predecessor, it only lowers its 
velocity before the predecessor escapes due to the kink- 
kink interaction and the successor becomes pinned at the 
boundary. 

So far we have determined the frustration strength 
X c (0) above which discommcnsurations enter the system 
in the absence of a driving force, / = 0. In the presence 
of a driving force /, we can roughly distinguish between 
the regime above equilibrium threshold, |x| > X<=(0), and 
below equilibrium threshold \\\ < Xc(0). Above equilib- 
rium threshold, \\\ > Xc(0), depending on the sign of x, 
kinks or anti-kinks are present in the system, since the 
entry barrier for discommensurations has vanished. Ne- 
glecting the effects of the Peierls-Nabarro barriers, the 
threshold force has basically vanished, 



/c«0. 



(26) 



For \x\ < Xc(0) there are no kinks present in the system 
at equilibrium due to the finite entry barrier. However, 
at a sufficiently large driving force, the entry barrier van- 
ishes, too. Then discommensurations enter the system at 
one boundary, freely flow through it and exit at the other 
boundary. In presence of a force the formation of kinks 
is similar to the kink-p enetration at equilibrium we dis- 
cussed in Section III C| . Of course, the extremal solutions 
which determine the energy barrier for kink formation 
are different. The stable and unstable solutions of the 
sine-X*ordon model in the presence of a driving force / 
readEI 



a$n + q 1 arcsin(///i) 



and 



(2n + l)ao 



i(//M), 



(27) 



(28) 



respectively. 

For frustrations close to the equilibrium threshold frus- 
tration, Xc(0) — |x| -C Xc(0), the dcpinning threshold 
force is small, f c <C fJ.. Let us study the depinning of an 
anti-kink at the left boundary x = for negative frustra- 
tions, x < 0, in the presence of a small force, < / <C fJ.. 
In the low-force limit, we can neglect the deformation of 
the anti-kink due to the force. The energy of the driven 
anti-kink relative to the uniform solution as a function of 
the anti-kink center x r is then 



AV(x c ) w — [u a JO; x c ) - u a , n (L; x c )] 
a 



(29) 



f 







j dx< 





dx 



AV s: 




FIG. 

A Vsan 



Potential energy of a small amplitude nucleus 



in units of V sa n = 4ka s /15aoh as a function of the 
position x s of the maximum amplitude. Shown are values 
close to the system boundary at x = for frustration param- 
eters x — — 2v3Xc(/)/5 (upper curve) and % = X<=(/) (lower 
curve), where the barrier vanishes. 



If we neglect the presence of the other boundary at x = L, 
the derivative of the potential with respect to x c is 



dx, 



■AV{x c ) 



K <9w a ,n(0; X c ) 



«0 



dx. 



X - 



duq, n (0;x c ) 
dx. 



Ua,n{0]X c ) 



(30) 



The anti-kink depins at the left boundary when the max- 
imum slope of the potential at x c = vanishes. This 
occurs at the threshold force 



fc 



4// 



1 



Xc(0) 



(31) 



which is easy to show realizing that u aj „(0;0) = ao/2, 
du a . n (0- 1 0)/dx c = Xc(0), and p = nk X 2 c \o)/(2ao). The 
same result is found for kink depinning at the other 
boundary, x = L, for positive frustration, \ > 0. 

At low frustration, |x| <^ Xc(0), the depinning thresh- 
old f c is close to /i. One thus has to consider the large 
force regime, [i — f <C where the lowest energy saddle- 
point solution of the sine-Gordon model u san ^ n (x; x s ) — 
u SiTl + Au(x]X s ) has a small amplitude and hence is 
called small amplitude nucleus (SAN). It can be calcu- 
lated by approximating the tilted cosine potential by a 
cubic parabola 



Am(i; x s ) — a s cosh 



2L 



(32) 



with center x s , amplitude a s = 3q 1 [2(1 — //^i)] 1 / 2 , and 
width I, = l d [2(l - f/v)]- 1/A . 

Let us now consider the penetration of an anti-kink 
for % < at the boundary x = 0. The energy difference 
AV S an = V[u sa n,n] - V[u Sjn ] as & function of x s < is 



AV san (x s ) 



15a l s 



3Aw(0, x s 



2a. 



1 



1 - 



_5 x Au(0,x s ) 

2^3 Xc(/) a s 
Au(0,x s ^ 3/2 ^ 



(33) 
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0.4 0.6 0.8 1 1.2 

X / z c (0) 

FIG. 6. Critical force f c as a function of the mismatch pa- 
rameter x- Shown are the numerical integration results (solid 
line), high force approximation f c — /t <C fi (dashed-dotted 
line) and the low force approximation f c <C /i (dashed line). 



where 



XcU) 



2^3 



2 1- 



3/4 



(34) 



As shown in Fig. |[ the SAN potential AV san (x s ) has 
a barrier (upper curve), which vanishes at x = ~Xc(f) 
(lower curve). This can be seen by analyzing the zeros of 
the derivative 



- — AV san (x s ) 
ax. 



Au(0,x s ) 



X 



ka s dAu{0, x s ) ( '< 
a ls dx s [3V3Xc{f) 

1/2) 



Au(0,x s ) 



(35) 



~Au~ 


2 


' Au~ 

a s _ 


4 


X 


a s _ 




~ 27 


lXc(f)\ 



One zero is given by x s =0, where the partial derivative 
dAu(0, x s )/dx s vanishes. The others can in principle be 
found by studying the term in the curly brackets which 
becomes zero if 



(36) 



However, to find the threshold condition, we do not need 
to calculate x s explicitely, it is sufficient to determine 
Au(0;x s ). Equation ([36j) has at most three roots, de- 
pending on the value of x 2 > 0. One of the roots is neg- 
ative, which is no solution, since < Au < a s [see Eq. 
([32|)]. For — Xc(f) < X < there are two positive roots, 
which indicate the existence of two extrema of A V san , a 
minimum and a maximum. Calculating the extrema of 
Eq. ( |3f^ ) with respect to Am, one sees that at x — ~Xc(f) 
the positive roots of Eq. ( p6| ) become degenerate, which 
means that the minimum and the maximum merge into 
a saddle point of AV san . Hence, the entry barrier of the 
SAN vanishes. At this value, Au = 2a s /3. Thus, two 
thirds of a SAN are localized at the left boundary, x = 0, 



but are unstable against small perturbations. Increasing 
either x or / depins the SAN which then evolves into a 
full anti-kink. Finally, from Au(0;a; s ) = 2a s /3 we find 
x s =l s ln(2- V3). 

To summarize, for low frustrations |x| *C Xc(0) and 
for frustrations close to the equilibrium threshold Xc(0) — 
|x| <C Xc(0) the threshold force f c is given by 



1 



Ixl 



Xc(0) 



4/3 



1 - 



Ixl 

Xc(0) 



0, 



Ixl « Xc(0) 

Xc(0) - Ixl « Xc(0) 
>Xc(0). 



(37) 



For completeness, we calculated the threshold force for 
arbitrary frustration by numerical integration. To calcu- 
late the static and dynamic solutions of Eq. (|l5|) , we use a 
standard numerical integration procedure. Starting with 
a flat initial configuration, u m = 0, we iterate 



u m (t + St) = u m {t) + St v m (t), 



with 



v m (t) = f + sin [27tm, 
- 2u m (t) + 5 rn 



(£)] + 2-Kg[u T 



l(*) 



(38) 



-i(*) 
(39) 



where length is measured in units of ao, time in units 
of to = Co^/Mj and force in units of fi. The length of 
the system is M = L/ao- Recalling that the vortex- 
vortex interaction energy falls off exponentially fast for 
distances between vortices larger than A, we take only 
the JV = [5 A/arj] next neighboring vortices into account 
in the sum over channel vortices. The channel has a 
length of L. At its ends, we apply boundary conditions 
taking into account the frustration, x- For a given force, 
Eq. (Bq) is iterated until a fairly steady state is reached, 
[v m (t + St) — v m (t)]/5t < 10~ 4 . Channel vortices are 
defined to be static, if v m < 10 -8 . The calculated u m (t) 
and v m (t) are recorded for several forces. In addition 
we can record the particle trajectories x m (t) to visualize 
the dynamical behavior of the channel vortices close to 
the depinning transition. The numerical results and the 
analytical limits for f c (x) are shown in Fig. ^|. 

The numerical integration also allows to determine the 
velocity averaged over time and space, 



L 

L ^ 



_ n jo 



dt 
f 



(40) 



as a function of the force / as shown in Fig. 0. 

Above Xc(0) where the entry barrier has vanished 
the topological defects move freely through the sample. 
The linear force- velocity characteristics resembles to the 
one of a single free particle with dissipative dynamics, 
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FIG. 7. Typical f — v curves computed for systems with 
L = lOOOao and A = do- The result for the unfrustrated com- 
mensurate case, x = (solid line), is identical to the curve 
of a single particle in a sinusoidal potential, v ~ (/ — fi) u 
with v = 1/2. Frustration lowers the threshold force / c , but 
does not alter the exponent v in the commensurate regime. 
However, in the incommensurate state, \\\ > Xc(0), the f — v 
curve becomes linear, v ~ /. Shown are the f — v characteris- 
tics for frustrated systems with \ = Xc(0) (dash-dot-dotted), 
0-7 Xc(0) (dash-dotted), and 0.5 Xc(0) (dashed). 

v = f/fj. The velocity of the entire chain is determined 
by the velocity of the defects that enter the system at the 
boundary. At |x| = Xc, where the entry barrier for de- 
fects has vanished, d%u(x = 0, L) = 0. This means that 
close to the maximum of the sinusoidal potential, where 
the chain bead spends most of the time, it effectively be- 
haves like a single particle which is driven by a force /. 
For |x| > Xc a similar argument can be given. At / = 
the effective force at the system boundary pushes defects 
into the channel until their density is so high, that their 
repulsion prevents new defects to flow in. Effectively, the 
chain bead at the boundary reaches an unstable equi- 
librium. Driving the system now with a non-zero force 
/ > results in the same motion for the bead at the 
boundary as for \x\ = Xc- 

For x < Xc(0) the entry barrier for the defects becomes 
relevant and the force-velocity characteristics shows the 
behavior typical for quasi-particles with a vanishing 
saddle-point barrier with v ~ [/ — fc(x)] 1 ^ 2 f° r / — * fc- 
Finally, in the absence of frustration the particles depin 
instantly and v = (/ c /r?)(2/// c - 2) 1 / 2 like a single par- 
ticle in a sinusoidal potential. 



E. Preliminary comparison with experimental data 

Let us now compare these theoretical findings with ex- 
perimental and numerical data.E3 Clearly, the general- 
ized Frenkel-Kontorova model predicts a nearly vanish- 
ing depinning force f c w in the defective state. In the 
commensurate state, the depinning threshold f c is finite, 
but decreases monotonously with increasing frustration 



and vanishes at the commensurate-incommensurate tran- 
sition, where the critical mismatch is reached, x — Xc(0), 
see Fig. ||. As a typical value for A Ian in the low field limit 
we consider A = 3ao. From Eq. (|13| ) we find g ~ 28 and 
hence Eq. @ yields Xc(0) ~ 0.06. Since X = (b /w) - 1, 
one expects a sharp spike in the threshold force around 
the magnetic induction with a vortex lattice plane dis- 
tance 6o = w - 

If we associate the maxima of the critical force to the 
spikes corresponding to a commensurate state and the 
minima to the incommensurate one, we would nonethe- 
less observe that some essential ingredient must be miss- 
ing in our model. Indeed, in Fig. 3 we see that instead 
of sharp spikes, the measurements rather show a con- 
tinuous modulation of the critical current as a function 
of the magnetic field, or equivalently, as a function of 
the frustration inside the channel. In addition, the data 
do not exhibit the vanishing pinning barriers which are 
expected for incommensurate (defective) structures with 
perfectly ordered channel edges. Both phenomena can- 
not be explained solely in terms of thermal or quantum 
fluctuations: They would lead to an effective reduction 
of the pinning barrier in both cases. This motivates us 
to study the influence of quenched disorder in the envi- 
ronmental vortex lattice. Actually, recent experiments 
have unvealed the essential role of disorder. It turns out 
that the picture of maximum critical force corresponding 
to an integer number of rows inside the channel is only 
valid for zero and weak disorder. If the disorder is strong, 
exactly the opposite will be realized, namely, the integer 
Bttmber of rows will correspond to the minima of f c (Ref . 

It is generally understood, that quenched disorder 
leads to pinning of vorticesa or topological defects of the 
vortex lattice, which would explain the increased criti- 
cal force in the defective state. It remains to understand 
how disorder leads to a reduction of pinning barriers in 
the commensurate regime. 



F. Channel with a distorted vortex lattice in the 
environment 

In a realistic sample the static nature of the channel 
environment is caused by some sort of pinning which may 
distort the vortex lattice. In the following we consider a 
device, where the vortices outside the channel are pinned 
by quenched disorder. The pinning in the channel envi- 
ronment is strong enough to guarantee that it remains 
pinned at all considered current densities whereas inside 
the channel pinning by quenched disorder is orders of 
magnitudes lower and can be entirely neglected. In a 
typical experiment the magnetic inductions are so large 
that the interaction between the vortices is much stronger 
than the pinning. Further, we assume that the vortex 
lattice in the channel environment is free of dislocations. 
Then it is natural to treat the vortex lattice in its elastic 







limit. In a weakly pinned vortex lattice, distortions of the 
order of the coherence length £ due to the disorder occur 
on average on a length l c in the direction of the mag- 
netic field and on a length R c transverse to it. Here, we 
consider two-dimensional collective pinning where l c 3> I 
such that the approximation of straight vortex lines re- 
mains valid. 

With quenched disorder in the environment of the 
channel, the interaction energy of a single vortex with 
the vortex row becomes 



E n {v) 



0. 



(41) 



where d m n are the displacements from the ordered posi- 
tions R m , n - It is convenient to rewrite the row potential 
in terms of the vortex density 



where 



E n (r) = J d 2 r'U{v-v') Pn {v'), 



p„(r') = ^2 S ( r ' - - d TO) „). 



We introduce a continuous displacement field 

,2 



d(r) 



4^ 2 Jbz tt 



(42) 



(43) 



(44) 



where BZ indicates that the integration is restricted to 
the first Brillouin zone. Note that d(R m-n ) = d mjIl . In 
order to derive a simplified one-dimensional model as in 
the ordered case, we make use of a few approximations. 
Since we consider the vortex lattice to be in the elastic 
limit, then 



|v-d| ~£/i? c « i. 



(45) 



Moreover, since the vortex potential falls off exponen- 
tially fast for |r — r'| > A, the channel environment is 
mainly probed within |r— r'| < A, where one can estimate 
|d(r')-d(r)| < \i)R c using Eq. @. For A£/i? c « a /2 
it is then reasonable to expand the displacements d(r') 
in the integral, 



d(r') = d(r) +0 



Rc 



(46) 



Here, we have introduced the coarse grained displace- 
ment field 

d(r) = {2ttU )- 1 \- 2 f d 2 r'U(r - r')d(r'), (47) 



which is smooth on the. scale A. Up to terms of the order 
0(\i/R c ), we obtainB 



E n {v)=E° n [v-d{v)] 



(48) 



for the collectively pinned vortex row potential. 

To calculate the effective channel potential for a chan- 
nel of width w ~ bo, we perform the summation over 
pinned vortex rows E n as in the ordered case. Further, 
since the influence of point-like disorder at the edge is 
assumed to be much weaker than the re-ordering due to 
the edge currents, we may well take d y (r) = 0. Then, 
the interaction of a single vortex in the channel with the 
disordered environment reads 



E dc (r) = A p (y) cos {q u [x - d x (r)]} 



(49) 



For the partial derivatives of the channel potential one 
finds 



d x E dc {r) = -2_^q u A u (y) sin {q v [x - d x (r)}} 
d y E dc (r) = ^A! v {y) cos {q v [x - d x (r)}} 

plus terms of order O (A£/i? c ). As in the perfectly or- 
dered case, we now consider the equation of motion. 
Substituting E oc by E dc in Eq. ( |ic| ) and following simi- 
lar arguments we derive an equation for the longitudinal 
motion. Introducing ip(x) = qd x (x, 0) we obtain a gener- 
alized phase-disordered Frcnkcl-Kontorova model, 



(50) 



n 1 J 



The corresponding energy functional in the continuum 
limit is then 



cos [qu — ip(x + u)]} — fu} 



(51) 



and the resulting equation of motion for the displacement 
fields u[x, t) reads 



rju = f — fj, sin [qu — tp(x + u)] + k 
+ kx[S(x -L)- S(x)}, 



d 2 u 
dx 2 



(52) 



where d x d x -terms that are of order 0(\£/R c ) have been 
neglected. 



G. Depinning in a channel with distorted 
environment 

Up to now, we have not specified the disorder displace- 
ment field d x . To gain some basic understanding, we now 
consider the effect of a local distortion on the depinning 
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properties. The field we choose is somewhat academic, 
but is convenient to understand the effect of a lattice 
distortion at the system boundary and in the bulk. The 
perturbation occurs around x v , 



threshold force for boundary depinning is again given by 
Eq. ([37]) with the modulus of the frustration replaced by 



W 

Ixl-y 



(59) 



d x (x,0) = W(x - x v ) 



9(x - x v ) 



1 



(53) 



Here, 9{x) is the Heaviside function and W > the dis- 
tortion parameter. It is convenient to introduce trans- 
formed displacement fields u = u — q^ 1 ^. Neglecting 
terms of the order (d x u) 2 <p' and ((p') 2 d x u, the relevant 
{t-dependent part of the energy functional reads 



To understand bulk depinning, we first restrict the con- 
siderations to x = and then discuss the behavior in the 
presence of frustration. 

For L/2 3> l s we can neglect the influence of the 
boundaries and treat them as if they were shifted to ±00. 
Then, the extremal threshold solution u v (x;x v ) can be 
constructed by joining two extremal solutions of Vsg- 
The matching condition is found from 



V[u] = V SG [u]+V x [u]+V v 



where 



V v [u] 



«o 



dx 



du dd x 
dx dx 



(54) 



(55) 



du(x — e) du(x + e) 



dx 



dx 



(60) 



such that the matching condition for the transformed 
field u v (x;x v ) reads 



In the following, we examine the effect of lattice dis- 
tortions at the boundaries, x v — 0, L, and in the bulk, 



diitpixif 



dx 



dx 



= W 



(61) 



= L/2. To gain a basic understanding of the depin- The mirror symme try requires d x u{x v ) = 0, hence 



ning process, we restrict the analysis to the large system 
limit, L/2 3> l s , where the system is so large that the 
depinning configurations at the weak spots x v — 0, L, 
and x v = L/2 do not interact with each other. 

For x v = 0, the lattice is distorted homogeneously in 
the entire sample and the contribution of the distortion 
to the energy functional yields 



kW 

V v [u] = ^-[u{L)-u{Q)]- 



(56) 



Physically, the distortion results in an additional frus- 
tration of the system, as can be immediately understood 
by comparing Eq. ( |5"6| ) with Eq. (|l7|). This means that 
for the threshold force, one can use the results that were 



d x u(x v ±e) = ^fW/2. At the threshold force f c stable so- 
lutions cease to exist. In fact, it can be shown that at f c 
the stable solution merges with an unstable one. This oc- 
curs, when the maximum of the tongue developing at x v 
reaches the maximum of the sinusoidal potential. ^From 
5V/5u — follows that at the potential maximum the 
extremal solution has to fulfill d 2 u(x v ) — 0, which holds 
if d x u{x v ± e) — ► for e — » 0. 

For weak distortions, \W/2\ <C I/C^n/s)) the threshold 
force for bulk depinning is close to fj, and the threshold 
solution u v = u s + Aii v can be found by merging two 
SAN solutions at x v ± x a , where x a = l s ln(2 + \/o), 



found in the absence of the distortion [see Eq. (37)], but 



Au lp (x;x ip ) 



Au(x; x v + x a ) x < x ¥ 

%/L {^X - X {p tJj (j J tX Xf 



(62) 



the frustration has to be replaced by an effective frustra- 
tion 



The maximum value of the tongue developing at x v is 
given by Au v {x v ;x v ) = 2a s /3. This implies 



X 



W 



(57) 



dAu v (x v — e;x v ) _ 2a s 



dx 



The result is thus a simple shift of the / c -curve. Sim- 
ilarly, for x v — L, the / — x-curve is shifted, 



3V3/ S 



W 

y 



(63) 



from which one obtains the bulk depinning threshold 
force in the presence of weak distortions, 



X -> X + 



w 

~2' 



(58) 



fc = M 



For x v = L/2, in addition to boundary depinning, bulk 
depinning at x v can occur. In large systems, L/2^> l s we 
can treat the effect of the distortions at the boundaries 
and in the bulk separately. At the boundaries, the thresh- 
old solution is then approximately given by the solution 
in absence of the defect, but with an increased effective 
frustration x+W/2 in the left half and a lowered effective 
frustration \ — W/2 in the right half of the system. The 



4/3 



(64) 



This formula becomes invalid for 2/{-K^fg) — \W\ <C 
2/{-K^/g). For strong distortions, \W\ < 2/(n^/g), the 
threshold configuration can be constructed by merging a 
kink and an anti-kink 



{^X 7 X^p J 



U'a 5 *^ tp ) ^ X^p 



(65) 
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from which one obtains the bulk depinning threshold 
force in the presence of strong distortions, 



fc 



4/x 



1 



*V9 



W 



(66) 



At W = 2/ (TTy/g) stable saddle-point solutions cease to 
exist for all / and disorder induced mobile kink-anti-kink 
pairs are spontaneously formed even at equilibrium. 

After having gained some understanding how pinning 
occurs at a weak spot in the bulk for \ — 0, let us now 
consider the frustrated case, |x| > 0, where bulk de- 
pinning competes with boundary depinning. Comparing 
the bulk depinning threshold f c (W) with the boundary 
threshold in presence of the defect, f c (\x\ — W/2), we find 
that for |x| < W the system depins in the bulk and for 
|x| > W at one of the boundaries. Note that to lowest 
order, we can apply these results to distortions of this 
kind that are not necessarily centered at x v — L/2, as 
long as l s <C x v <C L — l s holds. 



H. Depinning in a channel with randomly displaced 
edge vortices 

Opposed to the rather well behaved distortions of 
the previous paragraph, we now consider the effect of 
randomly displaced vortices in the channel edge. We 
mimic the disorder by uncorrelated relative displace- 
ments [i(Rm|i.„) — 4(Rm,n)]/"o- The latter are in- 
dependent identically box-distributed random numbers 
within the interval [- W/2, W/2], 

Pw{d x ) = ^ [8(d x + W/2) - 6(d x - W/2)] . (67) 

The width of the box distribution W parameterizes the 
disorder strength. Then, jl and d x (p are random func- 
tions, which are smooth on the scale A and bounded. 
Note that although d x (p is bounded, (p is unbounded. 
Thus, long range order is lost along the channel direction. 
On length scales much larger than ao, the displacement 
field d x (r) behaves like a random walk in ID and the 
phase-phase correlator scales linearly with the distance, 

The unfrustrated case has been discussed in Ref. 
Here we generalize the results for \\\ > 0. In weakly 
frustrated systems, |x| < W/2, nucleation still occurs at a 
weak spot in the bulk. Hence, the depinning threshold f c 
should be independent of the frustration \. This can be 
indeed observed, see Fig. |J the x~ /c-curve has a plateau 
around x ~ 0, before at larger values of x depinning takes 
place via the formation of defects at the boundary as in 
the ordered case. The threshold for boundary depinning 
is affected by lattic e dist ortions at the sample edge as 
discussed in Section III G resulting in an overall shift of 
the x ~ /c-curve to lower values of %. 

For intermediate frustration, <C x ^ Xc(0), a de- 
fect that entered the sample via the boundary cannot 
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FIG. 8. Numerically obtained minimum threshold force f c 
as a function of frustration \ for 100 channels with L = 100ao , 
A = ao, and W = 0.1. 
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FIG. 9. Static displacements u a (x) in a system of length 
L — 1000ao for x = — Xc(0) in the pinned regime, / < f c . (a) 
Pinned single interstitials (kinks), (b) Multiple pinned kinks 
at strong pinning sites, (c) Kink accumulation at the end 
of the system, x < L, where in addition to disorder the exit 
barrier is present. 



be pinned by disorder in the bulk. It moves freely to 
the other boundary, where it becomes pinned by the exit 
barrier until being released by the next defect that enters 
the channel and then travels freely to the exit. 

For frustrations above Xc(0) — W/2 bulk pinning be- 
comes possible since the boundary depinning threshold 
force becomes lower than the disorder induced pinning 
forces in the bulk. In Fig. || bulk pinning becomes rele- 
vant around x ~ 0.8x c (0). Indeed one can observe pin- 
ning of defects and bundles of defects in the numerical 
simulations in this regime. In Fig. |^ a static state in a 
system of length L = 1000ao for x — — Xc(0) just below 
the depinning threshold / < f c is shown. One clearly 
sees single pinned kinks and a few multiple-pinned kinks 
(bundles). At the exit of the system kinks accumulate 
due to both bulk pinning and due to pinning at the pres- 
ence of the exit barrier. When reaching the depinning 
threshold, a defect is formed at the entry on the channel 
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FIG. 10. Time evolution of the displacements u(x, t) in 
a system of length L = lOOOao for \ — ~ Xe(0) and / > f c . 
Strong pinning sites are indicated by bars parallel to the time 
axis. Here, interstitials (anti-kinks) travel to the right. At 
strong pinning sites they collide with pinned kinks. The for- 
merly pinned kinks are released while the previously moving 
ones become pinned. 
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FIG. 11. Disorder averaged f — v curves computed for 
systems with L — 100ao and A = ao- Shown are the / — v 
characteristics for systems with x — Xc(0) (dash-dotted), 
X = 0.5 ptc(0) (dashed), and \ = (solid). 



and travels until it collides with a defect that is already 
pinned at a strong pinning site. While the latter becomes 
released, the former gets pinned. This scenario repeats 
until a mobile defect has reached the channel exit, see 
Fig. 0. 

In contrast to the force- velocity curve of single systems, 
which behave as v ~ (/ — / C ) I/ with v = 1/2, the func- 
tional form of the disorder averaged force- velocity curves 
for finite systems all show upward curvature for / > f c , 
which corresponds to v > 1, see Fig. [H]. A crossover 
between exponents similar to the transition v = 1/2 — > 1 
occurring in the ordered model at Xc(0) can not be ob- 
served in the disordered case. 

Until now we considered the weak disorder limit in 
agreement with the assumptions made in order to develop 
the disordered model. If we increase the disorder param- 
eter W beyond the initially assumed limits, we can gain 
some insight into the depinning properties at large disor- 
der. In Fig. ^ minimum threshold forces f c as a function 
of disorder strength W for 100 channels with L — 100ao, 
and X = ao are shown for systems without frustration 
X = and for frustrated systems with x — Xc(0), where 
we expect the effect of the boundaries to become irrel- 
evant. At W — 0.5 the disorder is so strong, that a 
distinction between a commensurate and an incommen- 
surate system cannot be made. We speculate that this 
indicates a crossover to a depinning transition with true 
critical behavior as is reported for sandpile models or the 
Fukuyma-Lee-Rice model for charge density waves. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this work we have developed a model for 
artificial vortesJlow channels motivated by recent 
experiments.E§~Ej We have studied the depinning prop- 
erties of vortices in channels of finite length taking into 
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FIG. 12. Numerically obtained minimum threshold force 
f c as a function of disorder strength W for 100 channels with 
L — 100ao, and \ — ao- Upper curve: \ ~ 0, lower curve 
X = Xc(0). 

account inhomogeneities such as the sample boundaries 
and disorder. 

Throughout our analysis, we have neglected the influ- 
ence of thermal or quantum fluctuations on the depinning 
process. To see that this is well-justified with respect 
to the experiments that have been performed so far, let 
us estimate the magnitude of both types of fluctuations 
in conventional type-II superconductors which have been 
used to fabricate the channels. The typical pinning en- 
ergies of topological lattice defects are larger than the 
vortex self energy, which is about 10 eV, whereas the 
thermal energy is ~ 10~ 4 eU since the samples are cooled 
down to < 2K. Hence, only when the resulting energy 
barriers are reduced by a factor ~ 10 -5 one can expect 
to observe thermal creep of vortex lattice defects. Since 
the barrier energy scales as cx (l-j/jc) 3/2 , the width of 
the current density interval around j c at which thermal 
creep can be observed is Aj < 10 _3 j c . 

To estimate the relevance of quantum creep, we con- 
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sider the effective Euclidean action as a function of 
the current density, S e J^{j). At equilibrium, j = 0, 

S e J f U = 0) « & 2 /(e 2 p)Ud/jc) 1/2 , see Refl Here the 
coherence length is £ ~ 10 nm, the resistivity in the 
normal state is p ~ 10 _6 f2m, the depairing current is 

~ 10 9 A/m 2 , and the critical current is j c ~ 10 6 A/m 2 . 
We find S e J^ {j = Q)/h~ ~ 10 3 . For an overdamped vor- 
tex dynamics, which is characteristic for conventional 
superconductors, S^^(j)/h ~ 10 3 (1 — j/j c ) 3 ^ 4 - Quan- 
tum effects become irrelevant for S e J^ (j) ~ 7i, hence 
Aj ~ 10 j c . Thus, in artificial vortex-flow channels 
made of conventional superconductors thermal and quan- 
tum fluctuations only become relevant in an extremely 
narrow interval around the critical current which has not 
been resolved in the experiments discussed here. In a first 
step, we have therefore entirely neglected fluctuations 
and have studied a purely classical mechanical model at 
zero temperature. 

First, we have considered a perfectly ordered vortex 
lattice in the channel environment. Starting from a Lon- 
don description of straight vortices in a superconductor, 
we evaluate the periodic potential experienced by the vor- 
tices in the channel. The result, which is strictly valid 
only at inductions far below B C 2 is generalized to higher 
inductions by taking into account the vortex core and 
by introducing a field dependent coherence length and 
penetration depth. We have studied the overdamped dy- 
namics of a chain of interacting vortices in the effective 
channel potential. In our simplified description we have 
restricted the considerations to channels with a width of 
the order of the lattice spacing, w ~ bo. We have as- 
sumed that the resulting confinement perpendicular to 
the channel direction is strong such that the transverse 
motion of the vortices in the channel can be neglected. 
Further, since in all known channel experiments even for 
the lowest achievable magnetic inductions A > ao holds, 
to a very good accuracy only the first harmonic of the 
periodic potential is kept. These assumptions and re- 
strictions allow to determine the coefficients of a driven 
generalized Frenkel-Kontorova model. For A > ao, where 
the typical length Id of a topological defect in the chan- 
nel is much larger than the lattice spacing ao, but still 
much smaller than the system size l s , l 8 Id ao, the 
dynamics in the channel is conveniently described by the 
continuum limit of the Frenkel-Kontorova model known 
as the driven sine-Gordon model. 

After deriving the coefficients of the Frenkel- 
Kontorova model, we investigate the Commensurate- 
Incommensurate transition commonly known to occur 
in this model at thermal equilibrium. Since our focus 
is on the depinning problem in the absence of fluctu- 
ations, we modify the theory. We define a "mechan- 
ical" Commensurability-Incommensurability transition. 
It turns out, that the boundaries of the system play a 
crucial role if one supposes that the number of static vor- 
tices in the strong pinning environment is constant: the 
purely mechanical Commensurate-Incommensurate tran- 



sition occurs when the entry barrier for discommensura- 
tions at the boundary vanishes. 

This concept can be generalized to describe depinning 
at finite driving forces. The reason is that the entry bar- 
rier for discommensurations is reduced when a driving 
force is applied. The depinning then occurs at frustra- 
tions below the zero force threshold frustration, at which 
the mechanical Commensurate-Incommensurate transi- 
tion occurs. In this regime, the entry barrier is by far 
the largest barrier in the system. 

Above the zero-force threshold frustration, in the in- 
commensurate regime, the entry barrier has vanished and 
discommensurations enter the system until the mutual 
repulsion between the defects prevents new ones to flow 
in. The extremely small Peierls-Nabarro barrier which 
may arise due to the discreteness of the systemcJ is not 
taken into account. 

The depinning scenario depends on both, on the driv- 
ing force and on the frustration parameter X — { a ~ 
°o)/ a Oi which measures the mismatch between the lat- 
tice constant ao of the channel environment and the pre- 
ferred lattice spacing in the channel a. In the absence 
of frustration, % = 0, depinning occurs via a trivial ho- 
mogeneous solution when the barrier of the tilted wash- 
board potential vanishes and the threshold force is given 
by the amplitude of the sinusoidal channel pinning force, 
fc = l l - In the presence of frustration, depinning occurs 
via the formation of topological defects at one of the sam- 
ple boundaries. For weak frustration, \\\ "C Xc(0), we 
find f c = p{l - (1/2) [\/3 x/Xc(0)] 4/3 }, whereas close to 
the threshold frustration, Xc(0) — |x| "C Xc(0), we obtain 
fc = (4/x/?r)[l - x/Xc(0)] and for \ X \ > Xc(0) the depin- 
ning force vanishes, f c = 0. By performing a numerical 
integration of the equation of motion we have determined 
the values of / c (x) i n between these limits. 

Further, we have numerically calculated the force- 
velocity curves which correspond to the current-voltage 
characteristics of the sample. In the commensurate 
regime, \x\ < Xc(0) we find a characteristics typical 
to a saddle-point bifurcation scenario, v ~ (/ — fc) 1 / 2 , 
whereas in the incommensurate regime, where the com- 
mensurability gap has vanished, the response is linear, 
v ~ /. This behavior is typical for a system with open 
boundary conditions, where the density of topological de- 
fects is a function of both the frustration x and the driving 
force /. 

In systems with twisted periodic boundary conditions, 
where the density of topological defects is constant, the 
force-velocity characteristics is different in the incom- 
mensurate regimeEj. One observes a linear low-mobility 
regime for small driving forces, / -C p. The slope of the 
linear part is proportional to the density of topological 
defects in the system which is fixed by the magnitude of 
the twist at the ends of the system. The slope of the 
low-force regime is thus considerably smaller than in the 
high force regime / ^> p. For / > p the curves show a 
square-root dependence, v ~ (/ — p) 1 ^ 2 and only in the 
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high force limit they become linear again. The resulting 
force- velocity characteristics thus differ significantly from 
our results in the incommensurate phase, which is linear 
for all forces. 

We conclude that the boundary conditions strongly af- 
fect the force-velocity characteristics. The reason is that 
in frustrated systems the presence of (open) boundaries 
supports the formation of topological defects which lead 
to depinning. Further, the entry barriers at the bound- 
aries determine the rate at which defects enter the sam- 
ple, thus ruling the dynamic behavior entirely. The main 
problem in determining the behavior of vortex- flow chan- 
nels is thus to model realistic sample boundaries. Note 
that other boundary effects than considered here might 
modify the picture. For example, the vortex lattice of 
the channel environment may be distorted due to the 
presence of screening currents along the sample bound- 
ary causing a local variation of the frustration. Further, 
screening currents may lead to Bean-Livingston barriers 
for vortices which would increase the energy to form a 
defect at the sample edge. However, the conclusion that 
transport in artificial vortex-flow channels with a per- 
fectly ordered vortex lattice in the environment is ruled 
by the entry barriers at the sample boundary persists 
even if further boundary effects are taken into account. 

The picture that depinning occurs only via defect 
formation at the boundaries does not hold if inhomo- 
geneities are present in the bulk. For example, local dis- 
tortions of the vortex lattice in the channel environment 
caused by quenched disorder may generate weak spots 
that support the formation of vacancy-interstitial pairs 
at sufficiently large driving force. 

To investigate this issue, we have extended the model 
by accounting for small static displacements d(r) of the 
vortices in the channel environment. We assumed the 
channel environment to be in the elastic limit, |V • d| ~ 
£/i? c <C 1, where £ is the coherence length measur- 
ing the typical displacement within a distance given by 
the in-plane correlation length R c for lattice distortions. 
Further, we assumed that close to the channel edges 
A£/i? c <C ao/2 in order to obtain a local approxima- 
tion for the displacement fields of the static vortices close 
to the channel edge. Since £/i? c 1 this approxima- 
tion should be valid as long as the penetration depth is 
not orders of magnitudes larger than the typical vortex 
lattice constants. Further, we took only into account 
longitudinal displacements along the channel edge. For 
narrow channels, bo ~ w, we then obtained a general- 
ized amplitude and phase-disordered Frenkel-Kontorova 
model which in the continuum limit corresponds to a dis- 
ordered sine-Gordon model. Transverse displacements 
imply modifications of the disordered phase and ampli- 
tude of the sinusoidal pinning force and to additional 
pinning force terms. This issue and its consequences for 
depinning have not been considered here. 

In order to gain a basic understanding, we have first 
investigated depinning caused by a specific longitudinal 
vortex-lattice distortion field along the channel edges. 



Depending on their location, these distortions cause ad- 
ditional local frustration of the system modifying the 
threshold force for depinning. Distortions at the bound- 
ary of the sample affect the entry barrier for topological 
defects and cause shifts of the \ — f curve along the \~ 
axis. Local distortions in the bulk are shown to act as 
nucleation seeds reducing the threshold force. 

Finally, we have studied the effect of small disorder- 
induced displacements of vortices in the channel envi- 
ronment. We model disorder by uncorrelated relative 
displacements which are represented by random values 
that are independent identically distributed according to 
a box-distribution. In the absence of frustration local lat- 
tice distortions in the channel environment lead to an ef- 
fective channel potential with weak spots. At the weakest 
spots vacancy-interstitial pairs are formed when reaching 
the depinning threshold force. In the presence of frustra- 
tion, a crossover from bulk depinning to boundary depin- 
ning occurs when the entry barrier becomes smaller than 
the smallest bulk pinning barrier. Increasing the frus- 
tration the entry barrier is decreased until it becomes 
smaller than the typical bulk pinning force for kinks due 
to disorder. Applying a finite driving force that is large 
enough to overcome the entry barrier, but smaller than 
the bulk pinning force, topological defects enter the chan- 
nel and become pinned in the bulk. The finite depinning 
force in this regime is thus again determined by bulk 
properties. Increasing the driving force up to the thresh- 
old, topological defects which travel some distance un- 
til becoming trapped are successively introduced. Above 
the depinning threshold this leads to a jerky motion with 
successive depinning and pinning of topological defects. 

We obtained the force-velocity characteristics, the 
threshold force as a function of the frustration, and the 
threshold force as a function of the disorder for an ensem- 
ble of channels with randomly perturbed channel edges. 
Since the theoretical results are only valid for narrow 
channels with one mobile row, weak longitudinal disor- 
der, and since the squeezing effects of channel edge cur- 
rents were not taken into account, the interpretation of 
the available experimental data remains speculative. It 
would thus be interesting to extend the model to include 
these additional properties. 

Realizing that the Frenkel-Kontorova model and the 
sine-Gordon model are equivalent, the question arises, 
whether the disordered Frenkel-Kontorova model stud- 
ied herej-is-.related to the Fukuyama-Lee-Rice model 
(FLRM)E9c3 for charge density waves. This question 
is especially interesting with respect to the character- 
istics of the depinning transition. The FLRM and . sim- . 
plified versions haye-been studied both analyticalljStill 
and numerically^ L3 mostly in higher dimensions. The 
FLRM possesses a phase-disordered sinusoidal potential 
where the phases are chosen randomly from an interval 
[— 7r,7r]. This model shows critical behavior for d < 4. 
Approaching the threshold from below, the critical state 
is formed by the release of avalanches characterized by 
typical sizes that diverge with a power-law behavior. 
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Above threshold the motion is typically jerkyEHl and the 
velocity shows a power-law behavior v *~ (/ — f c ) v , where 
v depends on the dimensionality of the system. 

One is tempted to say that vortex flow channels pro- 
vide a physical realization of the one-dimensional FLRM. 
However, we find that the depinning process and the dy- 
namics above threshold strongly depend on the type of 
boundaries at the sample edges, the strength of the frus- 
tration, and the strength and type of disorder. In finite 
weakly disordered systems as studied in this work three 
depinning regimes can be identified. Increasing the frus- 
tration depinning first takes place via defect nucleation 
at weak spots in the bulk, then via defect nucleation at 
the boundary, and finally by releasing pinned pre-existing 
defects when the frustration is so strong that the entry 
barrier for defects has become irrelevant. 

This is indeed very different compared to the FLRM 
where system boundaries are not taken into account and 
the disorder is of a different type. If at all, characteris- 
tics of the FLRM like the avalanching below depinning 
threshold, the roughness of the critical state at threshold, 
and the jerky motion above the depinning threshold re- 
semble to what we observe in the incommensurate state 
in the presence of weak disorder. 

The main difference to the problem treated here is that 
the FLRM considers a system close to thermal equilib- 
rium whereas here we arc interested in the depinning 
behavior far from equilibrium. We identify typical con- 
figurations that act as sources which produce vacancy- 
interstitial pairs and thus lower the depinning threshold. 
These we call weak spots. We find also other configu- 
rations that pin the topological defects. The system is 
static if all topological defects that enter the system via 
the boundaries or are induced at weak spots are trapped 
by lattice distortions of the pinning type. Depinning 
takes place when the density of topological defects be- 
comes larger than the density of pinning sites. It is thus 
clear that an enhancement of the depinning threshold 
can only occur in systems with a considerable amount of 
pinning sites formed by lattice distortions in the channel. 
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